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Abstract 


Travel-times of acoustic signals were measured between a bottom-mounted source near 
Oahu and four bottom-mounted receivers located near Washington, Oregon, and Cali- 
fornia in 1988 and 1989. This paper discusses the observed tidal signals. At three out 
of four receivers, observed travel times at M2 and $2 periods agree with predictions 
from barotropic tide models to within +30 in phase and a factor of 1.6 in amplitude. 
The discrepancy at the fourth receiver can be removed by including predicted effects of 
phase-locked baroclinic tides generated by seamounts. 

Our estimates of barotropic M2 tidal dissipation by seamounts vary between 2 x 101% 
and 1 x 101? erg.s^!. The variation by two orders of magnitude is due to uncertainties in 
the numbers and sizes of seamounts. The larger dissipation (1 x 101? erg.s-!) is the same 
order as previous estimates and amounts to 4 96 of the total dissipation at M2. 
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Chapter 1 


Introduction 


We measure travel-times of acoustic signals between a bottom-mounted source near Oahu 
and four bottom-mounted receivers located near Washington, Oregon, and California 
in 1988 and 1989 (Figure 1-1). The power-spectra of basin-scale travel-times exhibit 
prominent variability at tidal periods [Spiesberger et al., 1989b]. There are two important 
reasons for further study of the tidal signals. First, numerical models for barotropic 
tides can be validated by using these new observations. Second, it is important to detide 
tomography data to better study other ocean fluctuations at similar periods [Spiesberger 
et al., 1989b]. 

The tides can be separated into barotropic and baroclinic modes, and their associated 
influences on travel-time can be separated into current and displacement effects. The tidal 
currents affect travel-time by alternating between alliance and opposition to the direction 
of sound propagation. The rise and fall of the water column affects the travel-time by 
causing periodic vertical iso-therm displacements and pressure fluctuations. Worcester 
et al. [1990] has shown that barotropic current effects can be quantified by using recip- 
rocal transmissions; one-way transmissions are used in this experiment, so we quantify 
barotropic effects by comparing measured data to barotropic tide model predictions. 

We find that acoustic travel-times computed from barotropic tidal models are similar 
to observations at three out of four receivers. The discrepancy at the fourth receiver 
can be understood by including the effects of phase-locked baroclinic tides generated by 


seamounts near the transmission path. Generation mechanisms are important both for 





A Seamount 





190 200 210 220 230 240775 


Figure 1-1: The acoustic tomography experiments of 1988 and 1989 consist of a source 
located near Kaneohe Bay, Oahu, and four bottom receivers (R1 to R4) whose approximate 


positions are shown. The distances between source and receivers are all about 4000 km. 





understanding the distribution of internal wave energy and the energy budget of tidal 
dissipation [Munk, 1966; Wunsch, 1975; Bell, 1975; Lambeck, 1980; Hendershott, 1981]. 
This report closely follows Headrick et al. [1990] and is organized as follows. Chapter 
2 discusses the experiment and the estimation of acoustic travel-time. Comparisons be- 
tween observations and predictions from barotropic and baroclinic tide models are made 
in Chapter 3. Chapter 4 predicts the global dissipation of tidal energy at M2 due to 


seamounts. Conclusions follow in Chapter 5. 
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Chapter 2 


Data Acquisition 


The travel-time differences of pulse-like signals are measured at 3 minute intervals in 
1988 and 2 minute intervals in 1989 between the source at Oahu and the four receivers 
(Figure 1-1). Travel-times are intermittently measured for a total of between four and 
eight months (Figure 2-1). The transmission distances are about 4000 km and the transit 
times are about 40 minutes. 

Spiesberger et al. [1989a,b] describes the maximum likelihood estimator of phase 
differences used to estimate differences in travel-time. The basic process is outlined in 
Appendix D. With this estimator, the change in acoustic phase at 133 Hz, 6¢, is related 


to the travel-time difference, 67, by, 


2764 
TU e : 
6 GOAL 0 


where w. = 27-133 (rad-s~!) is the center radian frequency of the transmitted signal, c is 
the radian frequency of the travel-time oscillation, and At is the geophysical time interval 


over which travel-time differences are measured (184 s in 1988, 122 s in 1989) [Spiesberger 
et al., 1989a]. 
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Figure 2-1: Shaded regions indicate the times when tomographic data were recorded at 


each of the receivers. 


12 





Chapter 3 


Data Analysis 


3.1 Barotropic Model for Acoustic travel-time 


The change in acoustic travel-time due to barotropic tides is given by, 
Aen) = = / ar hf OG as] " Y 
COMITE ር ES Co 

where the geophysical time at the receiver is t, the geodesic between the source and 
receiver is Ir, and the mean sound speed is c,. The barotropic current and surface 
displacement at location $ on the geodesic at time t are denoted by “(s,t — 5) and 
G(s,t — 9) respectively. The incremental on the geodesic is ds, and I(S) is the distance 
from the receiver. The constant is ky = 3.9 x 107? [Spiesberger et al., 1989b]. The effects 
of finite propagation time are taken into account by the correction term (4), The finite 
propagation delay amounts to a phase correction of order +10° at semidiurnal periods. 

The first term in (3.1) contains the effect of propagating sound with and against 
the tidal current. The second term, which amounts to a 1 or 2 ms correction, contains 
the combined effects of vertical iso-therm displacement and pressure fluctuations which 
accompany the rise and fall of the sea surface. For barotropic tides the travel-time change 
is about the same for different acoustic multipaths [Spiesberger et al., 1989a]. 

The elevation, ¢, is obtained from the Naval Surface Weapons Center (NSWC) tide 
model derived by Schwiderski [1978]. Unfortunately, except for M2, NSWC models of the 


tidal current field are unavailable. We therefore use a linearized version of Schwiderski’s 
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[1978] Laplace tidal equations (LTE) to predict tidal currents. The equations are, 


9u — 20 cos Ov + —L— |& (on ፦ ረ) — kgu, 


Ot Re sinô 
(3:2) 
at = -2N cos du 一 É [Zan — ረ) — kqv, 


where A and Ó are longitude and colatitude respectively, u and v are east and north 
current speeds respectively, 7 is Newton’s equilibrium tide (which is proportional to the 
astronomical tide generating potential), g is the earth’s gravitational attraction, Re is 
the radius of a spherical earth, 22 is the Earth’s angular rate of rotation, and kg is a 
damping coefficient. The constants o and f derived by Schwiderski [1978] correct for the 
solid-earth tide, ocean tidal loading of the solid earth, and the associated gravity potential 
perturbations of all tides and loadings. The solution to (3.2) for u and v is outlined in 
Appendix A. 

Equation (3.1) is decomposed into contributions from each tidal frequency, w,, of 
interest, and solutions are determined by collecting and integrating sinw;,t and coswy,t 


terms separately. The resulting travel-time change model is, 


K 


Al. Pr) = > ex(Tr) cos wyt + tell) sin wet, (33) 
k=l 


where e;(I'R) and f;(Tr) are the solutions to the separated integrals, and A is the total 


number of tidal frequencies calculated. 


3.2 Barotropic Model Comparisons with Observed Tidal 
Signals 


3.2.1 Tidal Signal Estimation 


Power-spectra of the differences in acoustic-phase exhibit prominent variability at tidal 
periods [Spiesberger et al., 1989a,b] (Figure 3-1). For each receiver, tidal signals are 
estimated by fitting the 8 most prominent tidal components to the data (Appendix A, 
Table A.1). Specifically, we determine a, and by; k = 1,2,...,8, in, 


8 
óp(tn, TR) = > al )cos art, + br Emeinirtn; n UNO... (3.4) 
k=1 5 
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i] en 


MT 


(cpd) 
Figure 3-1: Power spectral density estimate for the 133 Hz phase-differences at receiver 
R2. This periodogram uses a Bartlett data window with 42 days of data obtained in 1988. 
The tidal peaks at the 41, M2, and S2 frequencies have rms values which correspond to 


travel-time change amplitudes of 8.6, 23.5, and 12.3 milliseconds respectively. 


by linear least squares [Lawson and Hanson, 1974], where the phase-differences are mea- 
sured at times tn, and the total number of measurements equals N. Depending on the 
receiver, N equals 55,000 to 135,000, so the least squares problem is overdetermined. 
Each least-squares fit to phase-differences is compared to barotropic model travel-times 


to determine how much of the signal can be attributed to barotropic tidal effects. 


3.2.2 Comparison with Model 


Assuming all of the (3.4) signal is due to barotropic influences, each tidal component 
is proportional to the time derivative of the corresponding (3.3) model. Comparisons 
between observations and the barotropic model are illustrated by the phasors of Figure 


3-2. Each phasor magnitude corresponds to a data-to-model signal amplitude ratio, which 
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R1 R2 R3 R4 
Frequency Amp(ms) Amp(ms) Amp(ms) Amp(ms) 


M2 26 22 15 í 
52 11 T 4 3 
ለ1 10 T 1 3 


Table 3.1: Predicted acoustic signal amplitudes for barotropic tides. The model ampli- 
tudes, (e? + f?)!/?, are determined from the constants in (3.3) and have rms errors of 
about 2 ms. 


from (2.1) is, 
2r(a} + br)? 
wow, At(e? + f2)!/2’ 


and each phasor angle corresponds to a difference between data and model signal phases, 


msn cen 2 
tan e tan (3522 


which is given by, 


3.2.3 Analysis 


First, we discuss the comparison at two semidiurnal periods (M2 and $2). Predicted tidal 
amplitudes vary from 3 to 26 ms at our receivers (Table 3.1). At receivers R1, R2, and 
R3, observed amplitudes are 0.8 to 1.6 times predicted amplitudes (Figure 3-2). At R4 
the observed amplitudes are 2.4 and 3.4 times predictions. Predicted and observed phases 
are within +30° at R1, R2, and R3 but differ by up to 47° at R4. Data at R1, R2, and R3 
agree with predictions from the barotropic model to first order, but at R4, the barotropic 
model is inadequate. We return to this problem in the next section. 

We next discuss the comparisons for one semidiurnal period (X1). The models for 
barotropic currents at diurnal periods are inaccurate within a +10° band centered at 30° 
N (Appendix A). The transmissions to all receivers intersect this band, so comparisons 


between predictions and observations are problematic. Nonetheless, predictions at R1 
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Figure 3-2: Comparison of phase and amplitude of measured and predicted tidal signals 
at four receivers (R1 to R4) and three frequencies (M2, $2, KK1). Predicted tidal signals 
are obtained from the barotropic model given by (3.1). Each vector angle corresponds to 
a measurement phase that leads the predicted phase by the indicated amount. The unit 
circle in the upper left hand corner illustrates the vector angle convention. Vector magni- 
tude corresponds to the ratio of measured amplitude to predicted amplitude. Predicted 
magnitudes are listed in Table 3.1. An exact match between measurement and prediction 
is represented by a vector of angle zero and length one (tip just touching the unit circle at 
zero degrees). Solid arrows correspond to measurements of the combined 1988 and 1989 


record, and the dashed arrows correspond to separate 1988 and 1989 data measurements. 
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and R2 resemble observations at tidal period K1 (Figure 3-2). It may be that R1 and 
R2 tidal signals are similar to predictions because their transmission paths travel smaller 


distances within 30° + 10° N than paths to R3 and R4 (Figure 1). 


3.3 Baroclinic Tide Prediction 


3.3.1 Phase-Stable Internal Tides 


Apart from generation regions, baroclinic tides are not phase-locked to the ATGF, and 
their phases are incoherent at 100 km distance (approximately an internal tide wavelength) 
[Hendershott,1981]. In generation regions, baroclinic tides may be phase-locked to the 
ATGF [Munk et al., 1981]. Open ocean evidence of this phase-lock is seen in a phase- 
stable M2 baroclinic tide over a large guyot east of Hawaii [Noble et al., 1988]. The 
transmission path to R4 passes within a few hundred kilometers of several large seamounts, 
one of which is indicated in Figure 1-1. We hypothesize that the semidiurnal signals at 
R4 do not agree with predictions from the barotropic model because of the presence of 
phase-locked baroclinic tides near the seamounts. 

The transmission paths to R1, R2, and R3 may be near other large seamounts, but 
because there are no large discrepancies between the data and barotropic model predic- 
tions, it is unlikely that internal tides play a zero order role in determining tidal signals 
at these receivers. 

Appendix B summarizes a model for the generation of baroclinic tides from the inter- 
action between barotropic currents and seamounts. For baroclinic tides, acoustic travel- 
times are modulated by sound-speed fluctuations associated with vertical straining of 
iso-tachs [Spiesberger et al., 1989b]. The effect for the jth acoustic multipath, Pr(7), is, 

AT (STR) = af Gurt - DD (3.5) 
where the position along the multipath is 7, the vertical displacement of the internal tide 


is Gnt, and the vertical derivative of the potential sound-speed is der), 
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M2 92 
Amplitude(ms)  Phase(?) Amplitude(ms) Phase(?) 


Barotropic Model 7 113 3 324 
Baroclinic Model 20 117 5 265 
Combined Model 27 116 1. 286 
Measured Data 2፡1 1 DIES TOI 2009] 


Table 3.2: Comparison of predicted and observed travel-time changes at M2 and 52 for 
receiver R4. The uncertainties on the measured data are based on the differences between 
1988 and 1989 least-square fits. Phases are based on a sin wt time dependence with t — 0 
at 0000 Greenwich mean time on 1 January 1989. 


3.3.2 Influence on Travel-Time at Receiver R4 


Being the southern-most receiver, R4 has a more east-west oriented ray path than the other 
receivers. This orientation results in a weaker influence from the generally north-south 
oriented M2 and S2 barotropic tidal currents, so the effects of phase-locked baroclinic 
tides are more apparent. The predicted M2 baroclinic tide, generated from one of several 
seamounts near the transmission path to R4 (Figure 1-1), has an amplitude and wave- 
length of 6.2 m and 98 km respectively (Appendix B) and produces a 21 ms amplitude 
travel-time signal (Appendix C). 

The sum of the predicted M2 baroclinic tidal signals, generated from all prominent 
seamounts along the transmission path, is added to the barotropic prediction. The 
barotropic plus baroclinic prediction for R4 is almost identical to observations at M2 
and $2 (Figure 3-3 and Table 3.2). The prediction is improved at A2 but not at N2 
(Table 3.3). The frequency at N2 is close to M2 (Table A.1), and the observed amplitude 
at N2 is small. We suspect that the barotropic plus baroclinic prediction does not fit 
observations at N2, because strong baroclinic tides at Af2 may be smeared into N2 by 


dynamic processes [Hendershott, 1981]. 
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Combined 
Model 








Barotropic 
Model 





s, 
~, 
mn” 


Figure 3-3: The phasors on the left are identical to the combined 1988 and 1989 R4 record 
phasors shown in Figure 3-2. The phasors on the right present the same analysis using 


the combined barotropic and baroclinic models in Table 3.2. 


N2 K2 
Amplitude(ms)  Phase(?) Amplitude(ms) Phase(”) 
Barotropic Model 1 37 1 153 
Baroclinic Model 6 315 3 117 
Combined Model 6 325 4 126 
Measured Data cal 331.138 Sc 153 + 18 


Table 3.3: Same as Table 3.2, except comparisons are at N2 and K2 frequencies. 
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Chapter 4 


Tidal Dissipation 


The effect of the ocean tidal bulge on the Earth-Moon system is a well known phenomena 
that is discussed extensively by Lambeck [1980]. The friction-induced phase lag between 
C and 7 produces a torque in the Earth-Moon system that is responsible for the length- 
ening of the day and recession of the Moon at rates of about 1.8 ms-century”* and 3.7 


! respectively [Christodoulidis, et al. 1988]. The applied torque transfers an 


m-century” 
amount of orbital momentum to the Moon equal to the rotational momentum lost by the 
Earth (Lambeck, 1977). The net result of this momentum transfer is an overall reduction 
in system kinetic and potential energy. System energy is dissipated for the most part 
in oceanic tides, with some fraction (which has yet to be guantified) occurring through 


conversion to baroclinic tides [Lambeck, 1980]. In this chapter we estimate the barotropic- 


to-baroclinic conversion rate of the M2 tide at seamounts. 


4.1 Tidal Dissipation at M2 from Satellite Observations 


The rate of energy dissipation is, 


dE do d 
ーー(5.86ー- 9.58 )10erg - 71, (4.1) 


where HIP is the total tidal acceleration of the Earth and an is the acceleration of the 


Moon’s orbit in rad -s~? [Lambeck, 1977]. Recent estimates for M2 accelerations have 
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been obtained from the orbital perturbations of 17 satellites, and have values given by, 


OT - 445 x 10-2rad -s”?, CU 
Se = —9.73 x 107?frad s^ ?, 
[Christodoulidis, et al., 1988]. Substituting (4.2) into (4.1) yields 2.5 x 10!? erg-s~! at M2 


(about 75 % of all tidal dissipation). 


4.2 Tidal Dissipation by Seamounts 


We consider two models for the distributions of seamounts. In the first model the number 


of seamounts with a basal radius r > R (km) is given by, 
N(R) = 1.44 x 10° exp[—0.63R]. (4.3) 


Equation (4.3) is an extrapolation of Jordan et al.’s [1983] results, for the eastern Pacific, 
to the world ocean. In the second model, the number of seamounts is estimated from 
altimeter measurements by Seasat [Craig and Sandwell, 1988]. Seasat detects many large 
seamounts which are absent in the model by Jordan et al. [1983], but it does not detect 
small seamounts, as the altimeter profile sensitivity to seamounts falls off rapidly below a 
basal radius of 20 km [Craig et al., 1988]. These two seamount distributions are outlined 
in Table 4.1. 

The rate of energy conversion from barotropic to baroclinic tides at each seamount is 


approximately, 

E 三 和 和 zc Bis (4.4) 
where Az is the horizontal dimension of each seamount perpendicular to the direction of 
barotropic current, cg; is the group speed of the internal tide (Equation (B.22)), and F; is 
the average energy density of the internal tide (Equation (B.21)). The dimension, Az, is 
obtained from the seamount’s basal radius according to Table 4.2. If we take the barotropic 
current speed to equal 2 cm-s~! and the inertial frequency to equal 0.73 x 1074 s~! (30° 
N), then equations (4.4) and (B.23) yield a total energy conversion rate of 2 x 10!9 and 


1 x 108 erg-s~! for the exponential and Seasat distributions respectively. These estimates 


are proportional to Az and the squares of both seamount height and the barotropic 
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Number of Seamounts 


Basal Radius Exponential Seasat 
r(km) Model Model 
1-2 360, 000 0 
2-4 290, 000 0 
4 一 6 83, 000 0 
6-10 30, 000 200 

10 -- 15 2500 1200 
15 --20 110 6000 
20 -- 26 5 8800 
26 -- 32 0 6800 
32 -- 38 0 4800 
38 — 44 0 2800 
44 — 50 0 1600 


Table 4.1: Seamount distribution for the world’s oceans as a function of basal radius by 
Jordan et al. [1983] (exponential model) and Craig and Sandwell [1988] (Seasat model). 
The Seasat model is obtained by converting Craig and Sandwell’s [1988] histogram of 
peak-to-trough separation, p, to basal radius using r=(0.83)p in accordance with mea- 
sured seamount slopes. The numbers so obtained are multiplied by four to account for 
undetected seamounts [Craig et al., 1988]. Some of the discrepancy between models may 
be due to large seamounts in the western Pacific that were not represented in the Jordan 
et al. [1983] survey. 


Rectangular Seamount Dimensions 


Exponential Seasat 

Dimension Model Model 
Ar (0.62)r (0.62)r 
INS (0.21)r 2500 m 


Table 4.2: Seamount dimensions for the exponential and Seasat models (Table 4.1). The 
length and width of each seamount is Az, and the height is Az. Basal radius is denoted 
by r. The scale factor for Az, (0.62), is derived from a mean flatness of 0.31 reported for 
85 seamounts in the Pacific [Smith, 1988]. The Az scale factor for the exponential model 
is derived from the mean height-to-basal-radius ratio of 0.21 [Smith, 1988]. The height is 
set to 2500 m for the Seasat model. 
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current speed (Equation (B.23)). The exponential model yields insignificant dissipation 
because it contains relatively few wide seamounts (Table 4.1). The Seasat model yields a 
dissipation that is about 4 % of the total at M2. Munk [1966] and Bell [1975] obtained 
rates of 5.0 x 1018 and 2.5 x 1018 erg-s~! respectively for the conversion to baroclinic tides 
by bottom roughness. These estimates are the same order as we obtain using the Seasat 


distribution. 
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Chapter 5 


Conclusions 


We reach the following conclusions: 


1. Because barotropic signal predictions at Rl, R2, and R3 compare favorably with 
observations and because barotropic wave-lengths have basin-scales, the principle 
discrepancy between barotropic signal predictions and observations at R4 is probably 
not due to errors in Schwiderski’s [1978] model. Furthermore, because barotropic 
plus baroclinic signal predictions compare favorably with observations at R4, we 


believe the Schwiderski [1978] model is adequate for barotropic predictions. 


2. Seamounts may be responsible for a small (~ 5 %), but non-negligible, fraction of 
the dissipation of barotropic tides. The dissipation is sensitive to the heights, radii, 


and numbers of seamounts (all poorly known). 


3. The generation of baroclinic tides at seamounts can be accounted for with a two- 


layer model of the ocean (Appendix B). 


4. Barotropic tidal currents can be predicted from Schwiderski's [1978] tidal elevation 
predictions with our model (Appendix A), except within X10? of 30^ North and 


South (diurnal frequencies) and near 90? N (semidiurnal frequencies). 


5. Comparison of measured tidal signals (from basin-scale acoustic transmissions) with 
predicted tidal signals (from barotropic models) can be used to detect baroclinic 


tides which are phase-locked to the ATGF. 
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6. It is difficult to make a-priori predictions of tidal signals in basin-scale transmissions, 


because generation regions of baroclinic tides are poorly mapped. 
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Appendix A 


Obtaining Tidal Currents from 


Elevations 


We solve (3.2) for u and v by defining, 


u = bı coswt + b2 sinwt, 
: ' (ለ.1) 


ህ = ዐ1 ር08(፤ ተ ዐን 818 ሀ፤, 


and, 
eg ue de - = j 
Re | (om &c) = yı coswt + ysin wt, (A.2) 
-$ | (an _ ) = 21 coswt + z sin wt, 
where w is the tidal frequency of interest and a1, az, b1, b02, Y1, Y2, 21, and z2 are constants. 
The elevation, ¢, and Newton’s equilibrium tide, 7, are periodic with frequencies shown 
in Table A.1. The y's and z's are computed from (A.2) by expanding known functions, Ç 
and 7, in terms of sines and cosines. Substituting (A.1) and (A.2) in (3.2) leads to four 


linear equations, in the unknowns a1, a2, b1, and b2, which can be written in matrix form 


as, 


0 --90 cos O —w ka ay Ya 
—w E 0 20608 # ዐ 2 
_ IIA 
—22 cos 6 0 ka (の bi ፲1 
Ka の 20 cos O 0 ba 21 
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Semidiurnal Modes Diurnal Modes 


Frequency Frequency 
Component (cpd) Component (cpd) 
M2 1.93227 K1 1.00273 
52 2.00000 O1 0.92954 
N2 1.89600 Pi 0.99726 
K2 2.00547 Q1 0.89325 


Table A.1: Major Tidal Frequencies 


The solution to this matrix equation determines vu and v. 


We use the values, 

a 069. B= .90, (A.4) 
derived by Schwiderski [1978], and the value of the linear damping coefficient, kg, is 
selected by minimizing the phase difference between predictions of AT(t, M2, l'g3) using 
Schwiderski [1983] and damped LTE currents. The resulting best value is, 


ka 29.0 x 10-5s- 1. (A.5) 


The value of kg also controls the Q of the LTE model. We estimate Q by comparing 
the equations of motion to a forced linear-damped oscillator, 


ፅ*ፓ Ox 





387 +k 3 十 woz = Acoswt, (A.6) 
whose Q is defined as, 
WR 
“ሀ=” (ላ.7) 
where the resonant frequency is, 
k2 
wr = (wo — ie (A.8) 


[Marion, 1970]. By letting east be the z-direction, and substituting (A.2) into (3.2), we 
obtain the following z-direction momentum equation, 


"2 p 22 


37፤ + dT fv = yıcoswt + yasin wt, (A.9) 
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where f = 2Q cos @ is the inertial frequency. Since, 


|% - 11፡11. 
| fo| re (A.10) 
(7? + 32). = 4.4 ኗ 107". 
የ277 OU Mak, 
the dominant terms in (A.9) can be written 
E — fv - 0. (A.11) 
Substituting the solution to (A.11), 
y= aa (A.12) 
and the relation, 
U= WE (A.13) 
into (A9), we obtain the approximation, 
2 
22 + ka + wr = 91 COS の 7 + yasinwt. (A.14) 
Comparing (A.6,7,8) with (A.14) we see, 
(way o 1.41 x 10-4571 
2 (A.15) 


ka 0፡18. 5. 
where 1.41 x 1074s! is the frequency of the M2 tide. This agrees with estimates ranging 
from 5 to 25 provided by Hendershott [1981]. 

A comparison of the damped LTE model with the Schwiderski [1983] model for M2 
in the North Pacific ocean demonstrates good agreement (Figure A-1). Both models also 
generally agree with the data presented by Luyten and Stommel [1990] extracted from 
deep moored current meter records. By extension, we believe that the damped LTE 
model should be an accurate first order predictor of the $2, N2, and A2 semidiurnal 
currents. 

The damped LTE model does not perform as well as the Schwiderski [1983] M2 current 


model in the vicinity of amphidromic points. The finite difference approximations we use 
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Figure A-1: The M2 barotropic tidal current in the northeast Pacific Ocean at 0600 


200 


190 


Greenwich mean time on 1 January 1989 as determined from (a) the Schwiderski [1983] 


model (solid arrows) (b) linear-damped Laplace tidal equations (Equation (3.2); dashed 


arrows) (c) data extracted from deep moored current meter records (two arrows desig- 


nated by circles around their tails). The current meter analysis is taken from Luyten and 


Stommel [1990] 
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for the spatial derivatives in (3.2) are not accurate near these areas of rapidly changing 
phase. These inaccuracies are limited to a 100 km (the grid size of Schwiderski’s model) 
radius, and the net effect on predicted travel-time change is less than 5 % in amplitude 
and 3 degrees in phase for a 4000 km geodesic passing through an amphidrome. 

The damped LTE model has an additional difficulty for diurnal currents within a +10° 


band centered at the inertial latitude, 
0; = sin (wx/20). (A.16) 


For diurnal tides, w, ~ 1 cpd and 0+ = 30° N, and for semidiurnal tides, w, ~ 2 cpd and 
6; = 90° N. The difficulty can be seen by setting kg = 0 and solving (3.2) for u and v. 


The resulting equations, 





2.9 tw 0 
ም a|-20 A &(an-B6)| 
o R(4Q? cos? 0—4?) | 


(A.17) 





Fr an-Bl)-iw&(an- | 
በመ: 

have denominators equal to zero at the inertial latitude. These zeros are canceled by 
numerator zeros when u, v, and ( are determined exactly from the same equations of mo- 
tion, but a singularity occurs when the numerators are approximated by finite differences. 


The finite ky we use in (3.2) removes the singularity in (A.17), but the solution is still 


inaccurate near 30° N (Figure A-2). 
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Figure A-2: The A1 barotropic tidal current in the northeast Pacific Ocean at 0600 


190 


Greenwich mean time on 1 January 1989 as determined from the linear-damped Laplace 


tidal equations (Equation (A.2)). Note the irregular performance of the model in the 


N 


e 


vicinity of 30 
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Appendix B 


Guyot-Generated Internal Tides 


In order to determine the coupling between barotropic tidal currents and internal tides 
generated by guyots, we use the two-dimensional two-layer model illustrated by Figure B- 
1. This model is applicable to sea plateaus whose major horizontal extent is perpendicular 
to the direction of the barotropic current. We follow the procedure used in deriving 
the coupling between surface tides and coastal generated internal tides [Rattray, 1960]. 
Variables are a function only of the distance from the center of the plateau (y-direction) 
in this two-dimensional model. 


A rigid-lid approximation is used. The momentum equations are, 


00 + fu =0, (Bl 
Ov" 4 fu" 2 —9,0,C", (B.2) 
and, 
ðu! — fv' = 0, 
a (B.3) 


Zu” = fo" = 0, 
where u is the current parallel to the plateau, (’) terms correspond to upper layer and (”) 


to lower layer terms, f is the inertial frequency, and, 


gr = gAp/p, (B.4) 


is reduced gravity, where, 


App =p. (B.5) 
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Figure B-1: Definition Figure for a two-layer ocean model used to predict the coupling 
between a barotropic tidal current and the internal tide generated by a seamount whose 


extent perpendicular to the direction of the current is considered infinite. 
Subtracting (B.2) from (B.1), we find, 
dlv — v”) + f(u — u”) 2 g0,C", (B.6) 


and the continuity equations for each layer are, 


Q,(h^v') EC 


(B.7) 
の,( が の 7) = Scu 


Let vs = (a/d)expi(wt + 4), which is essentially constant for any distance y small com- 
pared to the wavelength of the surface tide. We have contributions to the currents in the 
two layers from both the barotropic current and the currents associated with the internal 


tide. Thus, 


の 5 us #u: Uu —u, tu! 


(B.8) 
v= vstu, u =u tiik, 
and from the rigid-lid approximation, 
ገው 
ae (B.9) 
E 
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We choose the following solutions which satisfy (B.3) and (B.7), 


y= 6 —l<y<é vet 
1 —Betkav Asin kıy Ber 
vs 1/d, 1/ሓ 1/d; 
v —(w[k;h'Be'* -(iw/kıh')Acoskıy —(w/kgh’) Be7**2¥ (B.10) 
v? (wk; T) Be'v (iw/kıhl)A cos kıy (w/kahl) Be**23 
us —tf/wd2 -if/wdı -if/wdz 
u (if[k;h')Be — -(ffkyh)Acosky (if /kgh') Be~**24 
u -(if/kahz)Be'®V (f/kaht)Acoskıy  —(if/kahz)Be"®v, 

with the common multiplier, 
aexp" 48). (B.11) 


To satisfy (B.6), 
= Sahn: aly: (B.12) 
In order to patch the solutions at the edges of the plateau, (" must be continuous at y = £ 
and y — —£. Thus, 
Beta! = Asin kyl. (B.13) 
The solutions (B.10) already satisfy the conservation of total mass flux h'v' -- h"v" — 
dv, + h'v: + h”v! across y = £ and y = 一 上 but conservation of p' mass flux requires that 


h'v' 2 h'(v, -- vj) be continuous also. To satisfy this requirement we let, 


Le Age a-r/2), 


| (8.11 
B- Bge- (07 7/2- 52). 
Substituting (B.14) into (B.13) leads to, 
Ao sin kıl = Bo. (B15) 
Using these relations and conservation of p’ mass flux we find, 
h' 5 h w ; . 
一 一 = Ao(sin a + icosa)sin kıl = — — — Ao(sin a + i cos a) cos kı £. (B.16) 
do ka d kh 
By equating the imaginary part of (B.16) we obtain, 
ki 
tana = y, n kit, (B.17) 
2 
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ly] <£ ly] > £ 


Amplitude Phase(°) Amplitude Phase(?) 
mt ዐ408]ኪ ደ17 $ — o t 90 aBo(|y|/v) 4 e -- 90 — kx(]y| — £) 
Us a/d; の a/do の 
v; awAocoskıy/kıh" ġ-a awBo/ksh' &- «+ 270 - ko(ly| — £) 
v; asAgcos kyy/ kh] | ó — o ተ 180 awBo/koh5 | $ — a 4- 90 — ks(|u| — £) 
Tin a f [cd $ 4- 270 a f [cda $ + 270 


u; afAocoskıy/kıh" | $ — a 4 210 afBo/k2h?  $— o4 180 — ka(|yl — £) 
ur af Aocoskiy/kih7 # - ወ ተ 90 af Bo/koh5 &—a— kollyl — £) 


Table B.1: Generalized two-layer internal wave solution for a barotropic north-south tidal 
current propagating over an undersea plateau of infinite east-west extent. 


and by eguating the real part, we can write, 


1 1 k,h'cosa 





DR ne a 8.1 
Vin በ ርር. (B16) 
which leads to, 
| eos 
Bo = (— — — ina. : 
0 (7 m っ sina ( B.19) 


The result of this solution is a standing internal wave over the extent of the plateau which 
patches to traveling internal waves at either end. The generalized solution is collected in 
Table B.1, and the particular solution for a guyot located at 30.3°N, 222.5°E is summarized 
in Table B.2. 

The cycle-average baroclinic energy density for this two-layer ocean is, 


— 1 


T 
Bi= | zPllCm)?g. + Wu)? + had? + hun)? + held, (8.20) 


where T is the period of the internal tide. Substituting the Table B.1 amplitudes and 


phases in (B.20) results in, 


Be Na 
== 2ዖ#("። 2 Bo) De ; (B.21) 
where the resulting energy density is in joules per square meter. This baroclinic energy 


propagates at a group speed that from (B.12) is given by, 


w? _ fe 
no DE (B.22) 
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M2 52 

Greenwich Greenwich 
Constituient Amplitude  Phase(?) Amplitude Phase(°) 
Vs 23 cms-1 68 0.7 cm-s~! 104 
vi 1.9 cm-s"! 346+ koy 0.6 cms- 22 4- kay 
vr 0.3 cms"! 166+ koy 0.1cm-s"! 202 +kay 
M. 12cms! 158 0.4 cm-s ^! 194 
ul. 1.0 cms"! 76+ ky 0.3 cms"! 112+ kay 
ur 0.2 cms"! 256 + koy 0.1 cm.s-1 292+ koy 


Table B.2: Two-layer M2 and 52 internal tide solutions for regions north of a 22 by 111 
km undersea plateau. Seamount is located at 30.3°N, 222.5°E (Figure 1-1), £ = 11 km, 
dj = 2200 m, dz = 4700 m, h^ 2 700 m, g- — .0059 m:s^?, k2(M2) = 6.41x107* m”! 
(wavelength is 98 km), and k2(S2) = 6.71x10~° m7! (wavelength is 94 km). 


The group speed, energy density, and horizontal extent of each seamount are combined 
in (4.4) to obtain the barotropic-to-baroclinic tidal energy conversion rate. Expansion of 


(4.4) in terms of fundamental parameters results in, 








2 /,2. ጋላጓ?"። 

一 = 2 (h' P! (dag, ) ^v? (1 一 7) (= "d | sin? a, (B.23) 

2 2 

where, 
/2 
a w- fr) 
t = t =- 一 一 Cle: ; 

ana (ans an o rh (B.24) 


The internal tide in Table B.2 has an energy density and group speed of about 1.5 x 10? 
erg - cm? and 160 cm-s”! respectively. The energy density near this seamount is about 
17 96 of that contained in the parent barotropic tide, and the baroclinic conversion rate 
is about 5.3 x 10!4 erg-s~! (53 Mw). 

The reduced gravity (B.4) used in the dispersion relation (B.12) is subject to the 
choices for the representative densities p' and o" in the upper and lower layers. Values of 


Gol order 10 2 ms” 


are typical for two-layer ocean models, but using g, = 107? m-s~? 
is not optimal for minimizing the mean square error between the combined model and the 
measured travel time signals at R4. The optimal value of reduced gravity chosen for this 


model, as illustrated by Figure B2, is about 5.9 x 1073 m-s~?. 
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Figure B-2: Magnitude of Receiver 4 barotropic plus baroclinic signal prediction errors. 


The errors are plotted as a function of reduced gravity. 
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Appendix C 


Baroclinic Tide Travel-Time 


Model 


Equation (3.5) predicts the travel-time modulation of a single multipath caused by one 
seamount. The multipaths oscillate about the sound axis, and a typical wavelength (upper 
plus lower loop) is 50 km. The potential sound speed gradient in the region of the Table 


B.2 plateau is approximately, 


0.13s^! above axis 


dz | 0.0s-! below axis, 


where positive z is upward. We assume the average ray path is above axis for 50 % of its 
travel adjacent to the 111 km length of the plateau, so the average change in travel-time 
at the M2 period is, 


(6.2 e" 2) +%cint m)(0.135-1)(111 km) 


. l 
= 51 tw(t—=-)+9¢, . | 
(2)(1479 m s-1)2 RL “ሠ... 


At) = 


where 04. is the phase of the M2 internal tide, and w is the M2 frequency. 

An evaluation of N, seamount-generated tides along the geodesic between source and 
receiver will produce N, baroclinic signals at each tidal frequency. The vector sum of 
baroclinic models plus the barotropic model results in the combined model, 


Ns 
ATi (tok) = ATpr(t,w;) + ATen tux), (53) 


n=l 
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where AT¿(n,t,w;) is the baroclinic signal prediction for the nth seamount at frequency 


Wp. 
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Appendix D 


Signal Processing 


D.1 The Transmitted Signal 


Pulse compression techniques described by Spiesberger et al. [1980] and Metzger [1983] are 
utilized to increase signal-to-noise ratio without degrading the phase difference resolution. 
The transmitted signal has a center frequency of 133 Hz and a bandwidth of about 17 Hz. 
The center frequency is phase modulated every 8 cycles using a 511 digit maximal shift 
register sequence having a period of 511(8/133) z 30.7 s [Spiesberger et al., 1989a]. The 


ideal transmitted waveform is, 
T 1l 
s(t) = so cos[27 fct — 3) - 3) (D.1) 


(Spiesberger et al., 1980), where fc is the center frequency, r, represents the shift register 


seguence of ones and zeros, and, 


8l 
一 一 )， 


133 (D:2) 


k(t) = 2 u(t — 
0 


where u(t) is the unit step function. 
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D.2  Post-Detection Processing 


D.2.1 Demodulation and Averaging 


The received signal, z(t), can be approximated by, 


Nmp 
ረ po basli = Ta); (D.3) 


n=1 
which is the weighted sum of Np multipaths each having gain br, and delay T, [Spiesberger 
et al., 1980]. Demodulation of the received signal is achieved by sampling at 4-133 Hz = 


532 Hz. Samples of an ideal single-path arrival are, 


amn] = cos(7-) cos[ (rita — 1/2) + sin( )sin[ (rut — 1/2)], (D.4) 
where, 
kin] = >= uln — 321], (D.5) 
1-0 


and u[n] is the unit step sequence. Because the sampling frequency is four times the 
carrier frequency, the even numbered samples are the cosine of the demodulated phase 
angle (weighted by +1) and the odd numbered samples are the sine of the demodulated 
phase angle (weighted by +1). The demodulated phase angles are low-pass filtered by 


block averaging to achieve a sequence of complex demodulates defined as, 
en = 27) fry y 0:1. .J = 1, (D.6) 


where the square root of minus one equals 2, 


2ሀ] = Do (-1)"*?2[2(m + 49)], (D.7) 


m=0 


and, 


3 
yli] = 9; (71)? z[(m 4 45) + 11) (D.8) 


m=0 


[Spiesberger et al., 1980]. The resulting sample rate is four complex demodulates per 
digit, so there are J = 4-511 = 2044 complex demodulates per sequence. The complex 
demodulates are averaged over 6 consecutive sequences in 1988 resulting in one averaged 
sequence record every 184 seconds. The demodulates are averaged over 4 sequences in 
1989. The signal-to-noise ratio increases by 10log,;,6 = 8 dB following the six-period 


average in 1988 and by 10log,,4 = 6 dB following the four-period average in 1989. 
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D.2.2 Replica Correlation 


One period of a 66.5 Hz sampling of (D.1) phase angles (expressed as complex demodu- 


lates) is contained in the following sequence of length J, 


T 1 T 1 
Sos cosl テ GAg ー 37] + ፤81812(፡:ሀ1 — 3); ld La (D.9) 
where, 
kb] = Yuli - 4l], (D.10) 
1-0 


and u[j] is the unit step sequence. The received complex demodulates (D.6) are circularly 
convolved with a sequence-reversed (matched filter) replica of the transmitted demodu- 


lates (D.9). The pulse-compressed output is given by, 


Ji 
dm] = = Sus ገ) በሙጭ01,..:,2..1. (1) ||) 
7=0 


where the periodic sequence, Sgljl, is given by, 


oo 


&sp]- 2, seD rJ], (D.12) 


[Oppenheim and Schafer, 1989], and s; is the complex conjugate of sy. The replica 
correlation (D.11) compresses 30.7 s of energy along each multipath into a single pulse 
of one digit length (~ 0.06 s) [Spiesberger et al., 1989a], whose amplitude and location 
depend on the gain, 6,, and delay, T,, in (D.3). The signal-to-noise ratio increases by 
1010፪16 511 = 27 dB following pulse compression, and a sliding average of 4 complex 
samples (separated by 0.015 s) increases the signal-to-noise ratio by another 6 dB at the 


cost of some resolution [Spiesberger et al., 1989a]. 


D.2.3 Phase Change Between Sequence Records 


The change in phase of the mth pulse-compressed demodulate from record r — 1 to r is 
given by the phase of, 
d[m,r — 1]d*[m, r]. 


The average change in phase from record r — 1 to r is given by the phase of 


M 
Aest = Y” d[m,r — 1)d"[m, r], (D.13) 


m=0 
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where there are M demodulates containing a signal [Spiesberger et al., 1989a]. The phase 
in (D.13) is referred to as the cross-correlated phase, because it is the phase of the complex 
cross correlation of adjacent pulse-compressed records at zero-lag. The cross-correlated 
phase can be converted to the average travel-time difference between records using (2.1). 
The precision of the resulting estimation of travel-time difference is about 135 js [Spies- 


berger et ዐ[., 108031. 
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